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We report updated results for Bk using HYP-smeared staggered valence quarks on MILC asq- 
tad lattices based on an analysis using SU(3) staggered chiral perturbation theory. The most 
important new feature of our data sample is the inclusion of a fourth ("ultrafine") lattice spac- 
ing. This improves the control over the continuum extrapolation and errors due our use of 
one-loop perturbative matching. We present a complete updated error budget, which leads to 
B K (NDR, n = 2 GeV) = 0.5309 ±0.005 1 ± 0.0424 and B K = Bjt(RGI) = 0.727 ± 0.07 ± 0.058. 
The results of the SU(3) analysis are inferior to those based on SU(2) staggered chiral pertur- 
bation theory, primarily because of the dependence on the Bayesian priors we use in the SU(3) 
fits. 
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1. Introduction 

This paper is the second in a series of four providing an update on our determination of Bk 
using improved staggered fermions. Here, we update the results obtained using fit functions derived 
from SU(3) staggered chiral perturbation theory (SChPT). In particular, we focus on the progress 
since last year's lattice proceedings [[]]]. Our present results are based on 10 ensembles of MILC 
asqtad lattices, whose properties are listed in Table [jj. As shown in the table, in the last year we 
have substantially increased the number of measurements on two ensembles, and added two new 
ensembles. In particular, the addition of the ultrafine ensemble Ul means that we now have 4 lattice 
spacings, as opposed to 3 last year. 

During the last year, we have also prepared a long article in which we give a detailed descrip- 
tion of all aspects of the calculation, and present final results based on 3 lattice spacings [g]. This 
reference contains all the details that we must necessarily skim over here due to space constraints. 
We also adopt the notation of Ref. [||] for parameters and fit types, and warn the reader that the 
labels for fits are different from those used in Ref. [|[]. 



a (fm) 


am(/am s 


geometry 


ID 


ens x meas 


B K (N-BB1) 


B K (N-BB2) 


0.12 


0.03/0.05 


20 3 x 64 


CI 


564 x 1 


0.555(12) 


0.564(17) 


0.12 


0.02/0.05 


20 3 x 64 


C2 


486 x 1 


0.538(12) 


0.535(17) 


0.12 


0.01/0.05 


20 3 x 64 


C3 


671 x 9 


0.562(6) 


0.592(14) 


0.12 


0.01/0.05 


28 3 x 64 


C3-2 


275 x 8 


0.575(6) 


0.595(13) 


0.12 


0.007/0.05 


20 3 x 64 


C4* 


651 x 10 


0.564(5) 


0.598(13) 


0.12 


0.005/0.05 


24 3 x 64 


C5 


509 x 1 


0.567(10) 


0.588(19) 


0.09 


0.0062/0.031 


28 3 x 96 


Fl 


995 x 1 


0.535(9) 


0.539(12) 


0.09 


0.0031/0.031 


40 3 x 96 


F2 # 


678 x 1 


0.540(8) 


0.545(13) 


0.06 


0.0036/0.018 


48 3 x 144 


SI* 


744 x2 


0.535(6) 


0.560(11) 


0.045 


0.0028/0.014 


64 3 x 192 


Ul* 


305 x 1 


0.540(6) 


0.547(8) 



Table 1: MILC asqtad ensembles used in the calculation. Ensembles marked with a * have improved 
statistics compared to last year, while those marked with a # are new. Results for Bk(h = 2 GeV) using both 
N-BB1 and N-BB2 fits are given. See text for discussion of these fits. 

2. SU(3) SChPT Analysis 

SU(3) SChPT was developed in Refs. f| and first applied to a calculation of B K in Ref. [§]. 
Since we use a mixed action, we need to generalize the SChPT calculation, and have done so in 
Ref. The result is that, for fixed a and sea-quark masses, the next-to-leading (NLO) order 
expression contains 14 low-energy coefficients (LECs). To obtain a good fit we also need to add 
a single analytic NNLO term, so that there are 15 LECs in all. Of these, 11 are due to lattice 
artifacts — either discretization errors or errors due to our truncation of the matching factors at one- 
loop order. The other 4 LECs remain in the continuum limit. 

On each ensemble we have 10 valence quark masses (running from ~ mf ,s down to ~ 
ra!( hys /10) and thus 55 different kaons. Nevertheless, a direct fit using all 15 parameters is not 
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stable, primarily because several fit functions are similar. To proceed, we reduce the number of 
parameters to 7 by removing 8 lattice artifact terms which have similar functional dependence to 
one of the 3 such terms that we keep. The details of this procedure are described in Ref. [Q]. This 
then allows stable fits. However, we find that, in some of these fits, the coefficients of the remaining 
lattice artifact terms come out larger than one would expect based on naive dimensional analysis. 
Thus as a second modification we constrain the size of the coefficients of these coefficients. We try 
different schemes for these constraints, as explained in Ref. 

We focus here on what we consider our most reliable SU(3)-based approach. This is based 
on a two step fitting procedure, in which we first fit only to degenerate points (m x = m y ), and then 
use the results of this fit as constraints on parameters for a fit to the entire data set. This approach 
makes sense because the the fitting form is much simpler for the degenerate kaons, containing only 
4 parameters (3 continuum and 1 lattice artifact). Fits to our 10 data points are stable and do not 
require Bayesian constraints, although we include such constraints for consistency, as we do need 
them in the second stage of the fitting. 

In slightly more detail, we fit to the form 

C = fctf (2.1) 

i=l 

where the functions F; depend on the pion and kaon masses and are given in Ref. F 4 is the 
lattice artifact term, and we constrain the LEC c\ by augmenting the % 2 : 

%Lg =X 2 + Zprior ; Zprior = (C4 - a 4 ) 2 / o% (2.2) 

We set «4 = (since we do not have prior knowledge of the sign of C4), and use 64 Aq CD (oAqcd) 
for the "D-Bl" fit or 64 ~ Aq Cd g^ for the "D-B2" fit. These two choices assume, respectively, that 
C4 is dominated by either discretization errors or truncation errors in matching. As noted above, 
these constraints have little impact on the degenerate fits. 

In the second stage, we extend the fit to the full data set, using the fit function 

^non-deg = ^ ^ 
i=\ 

where the first 4 terms are the same as in Eq. ( |Q| ) and the rest are given in Ref. [§]. Of the 3 new 
terms, only F5 survives in the continuum limit. We augment the x 2 with 

^prior ^prior (1) ^prior (2) (2-4) 
4 

Zprior(l) = ^( c i~ a i) l®i i #prior (2) = 2^ C jl®j 
«=1 ;=6,7 

Here a, ± a, are the results of either the D-Bl or D-B2 fit, which feeds in the information from the 
degenerate fits. The other priors are 

a („k ^2 



G7 



A^ CD (aA QCD ) 2 for N-BB1 fit 

A? 9 r -kt nnn r. ' \A&) 



a qcd« s 2 for N-BB2 fit 



A^ CD (aA QCD ) 2 forN-BBl fit ^ ?) 



A oa:>a? forN-BR2fit 



These constrain the lattice artifact terms and have a significant impact on the fits. 
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Figure 1: Bk (obtained using one-loop matching) versus Xp (squared mass of pion composed of valence d 
and d) for the C4 (left) and SI (right) ensembles. N-BB1 fits are shown. The red diamonds show the raw 
data, while the blue octagon shows the result obtained after extrapolation to to the physical quark masses 
with all taste-breaking lattice artifacts removed. 

3. Fitting and Results 

In Fig D, we show results of the N-BB1 fit on the C4 and SI ensembles. Compared to last 
year |JIJ], we have increased the statistics by factors of 10 (C4) and nearly 3 (SI). We also plot 
the data in a new way, showing it as a function of the squared mass of the pion composed of two 
light valence quarks. This displays the extrapolation to the physical kaon more clearly. Despite 
the reduction in error bars compared to last year (particularly in the left panel), the fit form gives a 
reasonable representation of the data. In Fig|2], we show a similar plot for the Fl and F2 ensembles, 
the latter being new. Again, the fits are reasonable. The quality of the N-BB2 fits is similar on all 
ensembles. Fits for the Ul ensemble are described in a companion proceedings 

Results for Bk from both types of fit are given in Table [|. We find that, on each ensemble, 
the two fits are consistent within 2a, but that a is relatively large, so the central values differ be as 
much as ~ 5%. This indicates a significant sensitivity to the size of the terms representing lattice 
artifacts, whose values cannot be pinned down by our fits alone. 

A striking result of the fits is that the coefficient C5 is very small. This multiplies 

F 5 ~( m f- m ff/{ m :? +m f), ( 3.i) 

which is the sole continuum term contributing only for non-degenerate kaons. The expectation is 
that C5 should be of 0(1), but it appears to be more than an order of magnitude smaller. This result 
implies that, in the continuum limit, we can almost determine Bk at the physical, non-degenerate 
point, using only degenerate kaons. This gives further a posteriori justification to our two-stage 
fitting procedure. 

A concern with these fits is poor convergence of SChPT. One can see from the figures that the 
result in the chiral limit (obtained by extrapolating the degenerate points to Xp = 0), which is the 
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LO term in ChPT, lies substantially below the final extrapolated result for Bk (given by the blue 
octagons). For more detailed discussion, see Ref. [0|. 



4. Continuum Extrapolation 

We use the results from the C3, Fl, SI and Ul ensembles to do the continuum extrapolation, 
since all have approximately the same values for mi and m s (and, in particular, the same ratio 
mi/m s ). We then extrapolate to the physical values of these masses, based on the dependence 
seen on the coarse and fine lattices. Note that the non-analytic part of this dependence has already 
been taken into account by setting mi = m^ hys and m s = m^ hys in the SChPT fit forms. In fact, the 
remaining dependence on mi is very weak, as can be seen from Table [l| 

The expected dependence on a is due to both discretization errors of the form a 2 a", with 
n = 0, 1, 2 . . . , and truncation errors starting at order a 2 , with a s evaluated at a scale ~ \/a [||]. We 
assume that the former dominate, with n = 0, and correct this assumption by adding in appropriate 
systematic errors. Thus we extrapolate using both linear and quadratic dependence on a 2 , as shown 
in Fig. I 

The data are consistent with both fit forms, although the quadratic fits are somewhat preferred. 
The parameters of the quadratic fit are, however, implausible. We expect the relative size of the 
quadratic and linear terms to be ~ (aA) 2 , which, with A = 500 MeV and for the coarse lattices is 
~ 0.09. Instead, in the quadratic fits the linear and quadratic terms are comparable on the coarse 
lattices. Thus we use the linear fits for our central values, and take the difference with the quadratic 
fit as the systematic error due to continuum extrapolation. 

The uuncation error we estimate separately by assuming the missing terms in the matching 
factor have size 1 xa s (l/fl) 2 , as explained further in Ref. [f2|]. 
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Figure 3: B K (NDR, \i=2 GeV) as a function of a 2 (in fm x 100) for the N-BB1 fit (left) and the N-BB2 fit 
(right), showing linear and quadratic fits to a = 0. 
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error ( 


%) memo 


status 


statistics 


1.0 


N-BB1 fit 


update 


matching factor 


4.4 


AB^ (SI) 


update 


discretization 


3.1 


diff. of linear and quadratic extrap 


update 


fitting (1) 


0.36 


diff. of N-BB1 and N-Bl (C3) 


® 


fitting (2) 


5.3 


diff. of N-BB1 and N-BB2 (C3) 




ami extrap 


1.0 


diff. of (C3) and linear extrap 




am s extrap 


0.5 


constant vs. linear extrap 


© 


finite volume 


2.3 


diff. of 20 3 (C3) and 28 3 (C3-2) 


@ 


scale r\ 


0.12 


uncertainty in r\ 


@ 



Table 2: Error budget for B K obtained using SU(3) SChPT fitting. 



5. Error Budget and Conclusions 

In Table g we collect our estimates of all sources of error, noting which results have been 
updated from our article We refer to that reference for details of the unchanged estimates. 

Compared to our result based on 3 lattice spacings the statistical error has decreased (from 
1.4%), as has the error due to the matching factor (from 5.5%). Both decreases are due to our 
addition of a fourth lattice spacing. On the other hand, our estimate of the discretization error has 
increased (from 2.2%). This is because we now use a different, more conservative method, namely 
the difference between linear and quadratic fits, rather than the difference between the result on 
ensemble S 1 and the continuum value. 

The largest error is now the "fitting (2)" error, which is our estimate of the uncertainty related to 
the different choices of priors in the Bayesian fits. We obtain this error from the difference between 
the results of the N-BB 1 and N-BB2 fits on the coarse lattices. We could use the difference between 
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these two fits after continuum extrapolation, which would more than halve the error, but we are not 
sufficiently confident in the continuum extrapolation to do so. In particular, the difference between 
these two fits remains substantial on the SI ensemble, whereas one would expect it to decrease 
compared to the C3 ensemble, since lattice artifacts are substantially smaller. 

Combining the systematic errors in quadrature, our current result for Bk using SU(3) SChPT 
fitting is 

5*(NDR,/i=2GeV) = 0.5309 ±0.0051 ±0.0424, 

B K = B K (RG\) =0.7270 ± 0.070 ±0.0580, (3 ' ' 

where the first error is statistical and the second is systematic. This updates our result Bk(2 GeV = 
0.524 ± 0.007 ± 0.044 given in Ref. [0. Compared with the SU(2) SChPT analysis the statisti- 
cal error is smaller here but the systematic error is significantly larger. Overall the SU(3) result has 
a larger error (8% vs. 5%). Given the less straightforward fitting, the concern with convergence of 
SU(3) SChPT, and the larger error, we use the SU(3) analysis as a cross-check on our result from 
the preferred SU(2) analysis. The results are consistent. 
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